Notes on the minimal longitudinal dc conductivity of perfect bilayer graphene 



Jozsef Cserti 1 

'Department of Physics of Complex Systems, Eotvos University 
H-1117 Budapest, Pdzmdny Peter setdny 1/A, Hungary 

We calculated the minimal longitudinal conductivity in prefect single and bilayer graphene by 
extending the two methods developed for Dirac fermion gas by A. W. W. Ludwig et al. in Phys. 
Rev. B 50, 7526 (1994). Using the Kubo formula which was originally applied for spintronic systems 
we obtain a™ n = (>^ 7r /2) e 2 /h while from the other formula used in the above mentioned work we 
find o™ n = (4«//vr) e 2 /h, where J = 1 for single layer and J = 2 for bilayer graphene. The two 
universal values are different although they are numerically close to each other. Our two results are 
in the same order of magnitude as that of experiments and for single layer case one of our result 
agrees many earlier theoretical predictions. However, for bilayer graphene only two studies are 
known with predictions for the minimal conductivity different from our calculated values. Similarly 
to the single layer case, the physical origin of the minimal conductivity in bilayer graphene is also 
rooted back to the intrinsic disorder induced by the Zitterbewegung which is related to the trembling 
motion of the electron. 

PACS numbers: 81.05.Uw, 73.23. Ad, 72.10.Bg, 73.43.Cd 



Unusual and remarkable transport properties of 
graphene (single or stacks of atomic layer of graphite) 
have been proved by recent experiments^. Besides 
the unconventional quantum Hall effect, many new phe- 
nomena have been predicted and studied in graphene 
such as the Klein paradoxon^, the specular Andreev 
reflection^, the Josephson effect, a new electric field 
effect, the photon-assisted electron transport 8 , compos- 
ite Dirac fermions^, quantum dotsiS, the n-p junctional, 
the fractional quantum Hall effect 12 , the spin-orbit gapii 

In this work we study the longitudinal conductivity 
of graphenes which, according to experiments, takes the 
minimum values of orders of e /h. This is an intrinsic 
property in the sense that it persists even in perfect (im- 
purity free) carbon honeycomb lattice. Such a peculiar 
behavior can be related to the excitation spectrum of sin- 
gle layer graphene described well by a Dirac like disper- 
sion relation (Dirac cones) of two-dimensional massless 
Dirac fermionsia. 

Much theoretical efforts have been devoted to ex- 
plain quantitatively the observed minimal longitudinal 
conductivity. However, at the theoretical side differ- 
ent predictions exist for the value of the minimal con- 
ductivity er™" even in prefect single layer graphene. 
Actually, long before the experimental evidence of the 
minimal conductivity, it has already been considered 
theoreticallyi£ii£ii! and Ludwig et al. found different val- 
ues using two different approaches^. Similarly, in many 
recent worka 19 i 20 ' 21 i 22 it has been found that o™ n = 
(4/7r)e 2 //i, while Ziegler predicted 2 ^ a™ in = ne 2 /h, 
Falkovsky and Varlamov 24 obtained <7™ n = (n/2) e 2 /h, 
and Nomura and MacDonald obtain from numerical 
calculations of the Kubo formula <r™" = (4/7r)e 2 /ft. 
for short range scattering case and o~™ n — 4e 2 //i 
for Coulomb scattering case 2 ^. Recently, Gusynin and 
Sharapov have derived analytical results for the ac and 
dc conductivity of Dirac fermions in graphene 2 ^ which 
could be used for studying the minimal conductivity. 



For bilayer graphene first studied experimentally 2 ! by 
Novoselov et al. and theoretically 28 by McCann and 
Fal'ko, the situation is not better regarding the mini- 
mal conductivity. Much less theoretical works considered 
the minimal conductivity in bilayer graphene. Recently, 
Koshino and Ando have investigated the transport in bi- 
layer graphene in self-consistent Born approximation 29 
and they found that in the strong-disorder regime er™ ra = 
(8/ir) e 2 /h, while in the weak-disorder regime <J™ n = 
(24/%) e 2 /h which is six times larger than in single 
layer graphene. Similarly, Katsnelson has also calcu- 
lated the minimal conductivity in bilayers using the Lan- 
dauer approach 29,3 - and he obtained a different value 
a™ n = 2e 2 /h. 

For massless Dirac fermion (corresponding to the 
single layer graphene) Ludwig et al. calculated the 
conductivity^ using two definitions for that, the one is 
the Kubo formula (c™™) and the other (a™ n ) is a def- 
inition often used in the sigma model literature 31 . The 
two definitions yield two different results (although nu- 
merically they are close to each other) for the longitudi- 
nal conductivity of perfect single layer graphene, namely, 
cr™ m = n vJ i s {Ti/8)e 2 /h and a™ n = n v n s (l/ir) e 2 /h, re- 
spectively. Here the factor n v = 2 and n s = 2 correspond 
to the two valleys and the electron spin, respectively (for 
details see references before). 

As can be seen from the above reviewed literature, 
several predictions of the exact value of the minimal 
conductivity even in perfect single and bilayer graphene 
exist although they are consistent at least in order of 
magnitude. In this work we extend Ludwig et al. ap- 
proach to calculate the longitudinal conductivity for bi- 
layer graphene. We find that, similarly to the case of 
single layer graphene, the two definitions of the conduc- 
tivity used in Ref. [l^ yield also different results for per- 
fect bilayer graphene. More interestingly, we find that 
the conductivity obtained from both approaches is dou- 
bled compared to that for single layer graphene. Thus, 
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our results show that obtaining the exact value of the 
minimal conductivity is a rather subtle task and it needs 
further investigation. 

Regarding the technical details of our calculations we 
note that when we used the Kubo formula, instead of 
following the approach of Ludwig et al. we applied an 
equivalent method developed by Bernevig^ for spintronic 
systems. In this method the conductivity is calculated 
in 'bubble' approximation using the one-particle non- 
interacting Green's function for finite temperature. This 
approach results in the same conductivity for single layer 
graphene as that obtained from Ludwig et al. method 
and, as can be seen below, it is more convenient to ex- 
tend to bilayer graphene. 

In order to treat the single and bilayer graphene si- 
multaneously we start from the Hamiltonian given in a 
unified form as 1,28 



Hi 





(p x +iPy) J 



(Px ~iPy) J 







(1) 



where J = 1 for single and J = 2 for bilayer graphene, 
and g is a constant depending on J. In what follows, it 
is more convenient to use the following equivalent form 
of Eq. @J 



Hj = fi(p)er, 



(2) 



where <x = (a x ,a y ) are the Pauli matrices representing 
the 'pseudo spin' acting on the spinor states with com- 
ponents corresponding to the wave function's amplitudes 
at the two non-equivalent basis atoms in the unit cell of 
the honeycomb lattice, fi(p) = gp J 
cp is the polar angle of p, i.e., p 
p = |p|. The eigenstates of the Hamiltonian Hj with 
plane wave solutions e lkr have eigenvalues 



" J (cos Jip, sin Jip) and 
p (cos ip, sin ip) with 



E ± {k) = ±Q(hk), 



(3) 



where ft = y£l 2 (hk) is the magnitude of the vector 

The Green's function of a Hamiltonian H is defined by 
G(z) = (z — H) and its position representation reads 
G(r,r',z) = (r|G(z)|r'). Using the Hamiltonian J5J one 
finds 



G(r,r',z) = (r|G(z)|r'} 

d 2 k ik (r-r') z + n(k)o- 



-n 2 (k) 



(4) 



where z is a complex number. 

To calculate the conductivity from the Kubo formula 
we apply the result derived first by Bernevig— for spin- 



tronic systems using the Green's function Q: 

e 2 . ha{Illf(cj)} 
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and n^^aj) is the retarded correlation function obtained 
by analytic continuation lo — > u> + ii] of function ILij(u>), 
np(E) — l/(eP( E ~^ + l) is the Fermi function with Fermi 
energy /i, the indices take the values 1, 2, 3, and summa- 
tion over repeated indices is assumed. Here A is the area 
of the two dimensional system, e a/ 3 7 is the Levi-Civita 
symbol, and the two energy bands E±(k) are given by 
Eq. ©. This result is valid for ballistic systems (note 
that this result is slightly rewritten to make it more trans- 
parent and useful for further calculations). 

We now evaluate the conductivity for undoped single 
and bilayer graphene at zero temperature using Eq. JSJ. 
The Fermi energy for undoped graphene is fi = 0, i.e., the 
negative energy band is fully occupied and the positive 
one is empty at zero temperature. Thus, in Eq. (|5b|) 
Uf(E + ) — riF(E-) = —1 for all occupied states with 
Fermi wave number k. For dc conductivity we need to 
find <j(u>) in the limit u — > 0. Therefore, the integra- 
tion over k can safely be extended to infinity since the 
contribution comes from the integrand close to k as 
(see the denominator of Using Eq. (|5bl) and 

— > A J d 2 k/(2ir) 2 (no contribution comes from the 
term contains the Levi-Civita symbol e a /3 7 ) we find 



H-M = lim, / ± 



2g 3 J 2 k 3J ~ 1 



n^o+ Jo 2tt (cj + i v ) 2 - Ag 2 k 2J 



(6) 



The imaginary part of the above integral can be evaluated 
analytically and it yields 



lm{U xx (u;)} = 



Thus, from i|5a|) the longitudinal dc conductivity (per 
valley per spin) for single (J = 1) and bilayer (J = 2) 
graphene becomes: 



8 h ' 



(8) 



Note that the conductivity obtained from the Kubo for- 
mula is universal depending only on the type of graphene 
but not the constant g. 

We now apply the other definition used by Ludwig 
et al. for calculating the conductivity (see Eq. (55) in 
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Ref.lH which reads as 



a xx =- lim i 7 2 /d 2 rr 2 Ti[G(0 > r > *T?)G(r,0,-itj)] i (9) 
h »,^o+ J 

where the Green's function is given by Eq. Q. The con- 
ductivity a xx can be rewritten asi&Si 

8e 2 d 
o xx = — lim rj 2 —r\ K(q,j]), (10a) 

K(q, 7] ) = ^Jd 2 r e lqr Tr [G(0, r, ir?)G(r, 0, -ty)] , (10b) 

and the trace is taken over the 'pseudo spin space'. 

Using the Green's function given by Eq. (0} and per- 
forming the trace over the 'pseudo spin space' in l|10b|) 
we find 



d 2 k 0+0_+t/ 2 



(lla) 



f2± = n(k±q/2). 



(lib) 



After a simple algebra the integrand / in Eq. (|lla(l can 
be calculated explicitly: 



I = 



(k 2 + q 2 /4) - (kq) 2 008(7$)+?? 



COS $ 



. 9 2 (fc 2 +^±k q y / +?? 2 , 

fc 2 - g 2 /4 
^/(fc2+qV4) 2 -(kq) 2 



-(12a) 

(12b) 
(12c) 



Here $ is the angle between vector k + q/2 and k — q/2. 
Expanding / for small q we obtain 



g 2 k 



2.1 



2 2 T.2J-2 



-q z g z Jk 



J (g 2 k 2J + 3?y 2 ) - [g 2 (J + 1) k 2J - (J - 1) V 2 ] cos (2y> - 2 7 ) 
4( 5 2 fc 2J + r; 2 ) 3 ' 



(13) 



where y> and 7 are the polar angle of k and q, i.e., k = 
fc(cos if, sin ip) and q = g(cos7, sin 7). As can be seen 
below 7 will be dropped out in the final result. 

Integrating the integrand / in i|13fl over k and taking 
the derivation with respect to q 2 we have the following 
simple result: 



(14) 



where 7 has been dropped out after the integration over 
ip. Similarly, the first term disappears in / since it is 
independent of q. Now using Eq. 110a|) the conductivity 
(per valley per spin) becomes 



J e 2 



(15) 



Discussion: We calculated the minimal conductivity 
for perfect single ( J = 1) and bilayer (J = 2) graphene 
using the Kubo formula and the definition 10. Taking 
into account the two valleys (n v — 2) and the electron 
spin (n s = 2) we find from Eqs. (JSJ and (|15l) that the 
minimal conductivities take the universal values: 



2 h 
4J£ 

7T h 



(16a) 
(16b) 



Note that these results are independent of the constant g 
appeared in the Hamiltonian . Just like in the case of 
single layer graphene we also find different values for the 
conductivity when the two methods in Ref.^is extended 
to bilayer graphene. Moreover, the minimal conductivity 
in perfect bilayer graphene calculated from Eqs. (J5J and 
© is doubled compared to that in single layer graphene. 
One can also see that for perfect single layer graphene 
our result l|16b|l agrees with that obtained in Refs. ll^l20l 

EM 

As pointed out by KatsnelsonSi^, the reason for the 
different results l|16a|) and l|16b|) is related to the 'trem- 
bling' or oscillatory motion of the center of a free wave 
packet called Zitterbewegung which has recently been re- 
considered in Ref. I33I The Zitterbewegung induces an in- 
trinsic disorder in perfect graphenes which has indirectly 
been confirmed by Tworzydlo et ali 2 ^ by calculating the 
Fano factor of the shot noise in single layer graphene, 
and they found that it has the same value (1/3) as in 
disordered metals. 

In our calculation the energy dispersion of bilayer 
graphenes was assumed to be purely parabolic in k (see 
Eq. (JTJ). Koshino and Ando studied the conductivity 
when the energy dispersion includes both k-linear and k- 
square terms 2 ^ which results in a trigonally warped Fermi 
surface 2 ^. To clarify the role of the k-linear term it would 
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be interesting to extend our work along this line. 

Similarly, extending the theoretical efforts to study 
multi-layer graphenes would be desirable in order to 
have better understanding of the unusual transport in 
graphenes. Recently, the electronic states in stacks of 
graphene layers have been studied^. However, it is not 
clear whether in the case of trilayers the Hamiltonian can 
be described simple by Eq. QJ with J = 3. 

The single and bilayer graphene are zero-gap semicon- 
ductors but this gap can open up due to impurities or the 
presence of interaction. Then, graphenes mapped to the 
massless Dirac fermion gas transform to random mass 
Dirac fermion systems. In this case, to calculate the con- 



ductivity it is useful to start from the definition © as it 
was demonstrated by Ludwig et al»i£. Recently, Gusynin 
et al. proposed theoretically that microwaves could be an 
experimental tool to understand the role of the impurities 
and the presence of interaction in graphene systems^. 

In summary, our work shows that the physical expla- 
nation of the existence of the minimal conductivity in 
graphene systems still remains a theoretical challenge in 
the future. 
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